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We introduce phase space concepts to describe quantum states in a disordered system. The merits 
of an inverse participation ratio defined on the basis of the Husimi function are demonstrated by a 
numerical study of the Anderson model in one, two, and three dimensions. Contrary to the inverse 
participation ratios in real and momentum space, the corresponding phase space quantity allows for 
a distinction between the ballistic, diffusive, and localized regimes on a unique footing and provides 
valuable insight into the structure of the eigenstates. 
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The behavior of a quantum particle in a disorder po- 
tential depends significantly on the disorder strength. If 
the mean free path exceeds the system size, one may 
think of plane waves characterized by a fixed momentum 
which are scattered by the weak random potential. On 
the other hand, sufficiently strong disorder leads to ex- 
ponentially localized states in real space. Then a wide 
range of momenta is needed to construct these states. A 
description for arbitrary disorder strength thus requires 
to adequately take into account real space as well as mo- 
mentum space properties. In this paper, we will therefore 
adopt a phase space approach. 

Signatures of the different regimes can already be 
found in the energy spectrum. While for weak and strong 
disorder energy levels can almost be degenerate, level re- 
pulsion occurs in an intermediate regime. Making use of 
random matrix theory the chaotic nature can be verified 
and related to diffusive motion. However, such consid- 
erations cover only statistical properties of the spectrum 
and do not give information about the structure of indi- 
vidual states. 

A popular way to investigate the properties of single 
states is the calculation of their inverse participation ra- 
tio. In real space this quantity has frequently been em- 
ployed ||l[ to measure the size of the localization domain 
of quantum states in the localized regime and to char- 
acterize the Anderson transition |^. However, the 
real space inverse participation ratio is not very sensi- 
tive to changes in extended wave functions when going 
from the ballistic to the diffusive regime. In order to ob- 
tain a meaningful description of the structure of quantum 
states for all regimes on a unique footing, we generalize 
the concept of the inverse participation ratio to phase 
space. 

The Anderson model of disordered solids has been 
the subject of extensive investigations over the last 
decades A numerical study of this model in one, 
two, and three dimensions will demonstrate the virtues 
of our approach. In particular, we will be able to identify 
ballistic, diffusive, and localized regimes via properties of 



the eigenstates. In one dimension the results will signifi- 
cantly differ from those in higher dimensions, which can 
be attributed to the absence of a diffusive regime. 

We start by introducing the relevant phase space con- 
cepts. At this point there is no need to specify the details 
of the disordered system except that we will consider a 
d-dimensional lattice model with lattice constant a and 
length L in each direction. In order to keep the notation 
simple, we will give the formulae for the case of one di- 
mension, which can be generalized to higher dimensions 
in a straightforward manner. 

A positive definite density in phase space is given by 
the Husimi function |^ or Q function Q 



g{xo,kQ) = |(a;o, fcolV")!^ 



(1) 



Here, the state {ip) is projected onto a minimal uncer- 
tainty state |a;o, fco) centered around position xq and mo- 
mentum ko. In position representation, the latter as- 
sumes a Gaussian form 



(a:|a;o, fco) = 
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where the value of the variance is yet undetermined. 
The Husimi function is normalized, J(dxdk/2Tr)g(x, fc) = 
1, and, for real wave functions ip{x), obeys the symmetry 
g{xo, fco) = gixo, -fco). 

The variance of the Gaussian (^) determines the 
relative importance of real and momentum space. In the 
following, we choose <t^ — La/An leading to equal widths 
of the Gaussian in x- and fc-direction. In order to obtain 
sufficient resolution, one has to ensure that a L which 
limits the possible system sizes from below. Then, the 
effect of neglecting the tails of the Gaussian in finite size 
systems in presence of periodic boundary condition will 
also be small. Examples of Husimi functions for one- 
dimensional disordered systems are presented in Fig. 
which will be discussed in detail below. 

In higher dimensions, Husimi functions can no longer 
be visualized easily. A more global description of the 
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FIG. 1: The Husimi function of a state at tlie band center is shown for a one-dimensional Anderson model of size L = 128 with 
disorder strengths W — 0.1, 1, 2.5, 10, and 25 increasing from left to right. Exploiting the symmetry with respect to fc = the 
Husimi function is plotted on a linear and logarithmic scale in the upper and lower half, respectively. The gray values from 
white to black for increasing Husimi function have been normalized separately in each panel. The width a of the Gaussian (^) 
is indicated by the black square in panel (e). 



phase space properties of the state is therefore necessary 
and often sufficient. Based on the Husimi function the 
so-called Wehrl entropy |l^ is defined, which rep- 

resents a measure of the phase space occupation. It was 
shown for the driven rotor that the Wehrl entropy of in- 
dividual quantum states is connected to the energy level 
statistics A very similar system, the kicked rotor, 

can be mapped onto the Anderson model suggest- 
ing that the Wehrl entropy is a useful quantity for the 
characterization of the eigenstates of the Anderson model 

For practical purposes it is more convenient to intro- 
duce the inverse participation ratio in phase space 

-^[Qix,k)]\ (3) 



P = 



which corresponds to a linearization of the Wehrl entropy 
[Q . This inverse participation ratio can be compared di- 
rectly to the corresponding quantities Px — J dx\tlj{x)\'^ 
and Pk = J dk\^{k)\'^ in real and momentum space, re- 
spectively. Px is particularly well studied as it is related 
to the probability for a diffusing particle to return to its 
original position in the long-time limit [ p^ . 

Furthermore, P can be evaluated without recourse to 
the 2d-dimensional Husimi function. In fact, only the 
wave function ipi^) is required since we can recast in 
the form 



P = 



1 



80Fct 



du 



dvi^l ^{u~v)j + ) (4) 
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Its nondiagonal character provides the information on 
momentum. It is only by means of (Q) that one succeeds 
in determining the inverse participation ratio for three- 
dimensional systems. 

In the following, we specifically consider the Anderson 
model for non-interacting electrons on a lattice with pe- 
riodic boundary conditions, where disorder is modeled by 



a random on-site potential. In d = 1, the Hamiltonian 
reads 

H = -t^{\n){n+l\ + \n+l){n\) + W^Vn\n){n\ (5) 



with Wannier states |n) localized at sites n — 1, . . . ,L. 
All lengths are measured in units of the lattice constant 
a = 1. The hopping matrix element t — 1 between neigh- 
boring sites defines the energy scale. The on-site ener- 
gies Vn are drawn independently from a box distribution 
on the interval [—1/2; 1/2] and W denotes the disorder 
strength. 

Husimi functions for a state at the band center are 
presented in Fig. |^ for increasing disorder strength W 
and a randomly selected disorder realization w„ . Making 
use of the symmetry with respect to fc = 0, the Husimi 
function is plotted on a linear and logarithmic gray scale 
in the upper and lower half, respectively. White points in 
the lower half may be related to the zeros of the Husimi 
function [p^ . 

For W = 0.1 (Fig. |^a), the disorder represents only a 
small perturbation and the Husimi function is thus still 
close to that of plane waves. Except for the states at the 
band edges, one finds two stripes well localized at the 
corresponding fc-values, which are extended over the full 
real space. The width of the stripes is induced by the 
projection onto the Gaussian (||). In Fig. |l]e the opposite 
limit is depicted. At = 25, the hopping is a small 
perturbation, t ^ W, and the state is localized in real 
space. 

In Fig. I^d, where W — 10, the influence of the nearest 
neighbor hopping becomes relevant and tends to extend 
the wave function in real space over several sites. Since 
the coupled states have to remain orthogonal, they sep- 
arate in momentum space. This leads to a contraction 
of the Husimi function in k direction, which is more im- 
portant than the spreading over a few sites. As a conse- 
quence, the phase space properties for strong disorder are 
dominated by the behavior in momentum space and the 
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FIG. 2: Distributions of the logarithms of the inverse partic- 
ipation ratio for 50 disorder realizations in (a) real space, (b) 
phase space, and (c) momentum space as a function of the 
disorder strength W for a one-dimensional system of length 
L = 2048. 

inverse participation ratio in phase space increases with 
decreasing disorder. This behavior will also be evident 
from Figs. || and |3| below. 

Similarly, at weak disorder, the transition from Fig. |l|a 
to Fig. I^b, i.e. to W = 1, can be understood in terms of a 
coupling between different plane waves. Here, however, 
the coupling is not restricted to neighboring k values, 
but is governed by the energy difference of the respective 
states. For one-dimensional systems, the contraction in 
real space dominates the spreading in momentum space, 
again leading to an increase of the inverse participation 
ratio in phase space (cf. Fig. ||b). 

As a consequence of the behavior for weak and strong 
disorder, one expects a maximum for the inverse partic- 
ipation ratio at intermediate disorder strength. Indeed, 
for W — 2.5, the state shown in Fig. IJc displays strong 
localization in phase space. 

The situation just described is generic for one- 
dimensional systems. This can be seen from the dis- 
tributions of the inverse participation ratio depicted in 
Fig. Ij. The distributions of the logarithms of P as well 
as Px (real space) and Pu (momentum space) have been 
obtained by diagonalizing Eq. (^ for 50 different disorder 
realizations w„ for each disorder strength W and taking 
L/2 states around the band center into account. 

In the limit of very strong disorder, the states are lo- 
calized on a single site in real space and uniformly dis- 
tributed over momentum space. This leads to the limit- 
ing values Pxipo) = 1 and Pk{oo) = L~'^. In phase space 
one has to account for the finite width of the Husimi 
function and thus finds P{oo) = L^"^/"^. For d = 1, these 
values can be checked against the data shown in Fig. |^. 
Starting from this limit, with decreasing disorder two en- 
ergetically almost degenerate states become coupled via 
the finite hopping matrix element t. For these states Px 
is reduced to 1/2 while Pk is enhanced by a factor of 
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FIG. 3: Distributions of the logarithms of the inverse partic- 
ipation ratio in phase space (top) and r] defined in Eq. (^) 
(bottom) are shown as a function of disorder strength W for 
(a) d = 1, L = 128, (b) d = 2, L = 64, and (c) d = 3, L = 20. 
The distributions are based on 250, 20, and 20 disorder real- 
izations, respectively. 

3/2. As already discussed above, it is the latter which 
dominates the behavior in phase space. 

In the opposite limit W ^ 0, the real wave functions in 
d = 1 contain equal contributions from degenerate plane 
waves of momenta k and ~k. This implies Pfc(O) = 1/2 
and Px{^) = 2/3L. In phase space, the finite width of 
the Husimi function leads to P{0) = 1/2L^/^. In higher 
dimensions, however, degeneracies occur and render the 
behavior for W ^ more complex. Nevertheless, as 
a function of system size, Px and P scale as L~'^ and 
j^~d/2^ respectively. A more detailed discussion will be 
given elsewhere |l7| . 

The upper part of Fig. ^ shows the distributions of 
the logarithm of P for dimensions d — 1,2, and 3 and 
system sizes L = 128, 64, and 20, respectively. For better 
comparison the data have been scaled with the length 
dependence L^"*/^, being valid in the limits W ^ and 
oo. 

The most striking difference between the one- 
dimensional case and higher dimensions consists in the 
behavior of the inverse participation ratio for weak dis- 
order. For d = 1, the average inverse participation ra- 
tio increases with disorder strength and eventually goes 
through a maximum. The overall behavior can be un- 
derstood as a crossover from a regime dominated by real 
space properties to one dominated by momentum space. 

In contrast, in d > 2 the average inverse participation 
ratio initially decreases. This implies a spreading of the 
Husimi function beyond the broadened stripes present for 
weak disorder (cf. Fig. ^). On the other hand, the be- 
havior at strong disorder is governed by the same mech- 
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anism as for the one-dimensional case, which impHes an 
increase of the inverse participation ratio with decreasing 
disorder. The two regimes are joined by a short interval 
of disorder strengths where the Husimi function contracts 
strongly as disorder is increased. 

The distributions of P depicted in Fig. || behave sim- 
ilarly for d = 2 and 3. However, the pertinent scaling 
argument for Anderson localization |^ predicts a phase 
transition only in dimensions higher than two. Unfor- 
tunately, in d = 3 numerical constraints prevent us from 
performing a finite size scaling, which would allow to dis- 
tinguish the form of the jumps of P in d = 2 and 3 in the 
thermodynamic limit L — > oo. 

The main difference between the behavior in d = 1 
and d > 2, the appearance of a minimum in the inverse 
participation ratio P, can be associated with the exis- 
tence of a diffusive regime in d > 2. While in d = 3 
a diffusive regime appears even in the thermodynamic 
limit, in d = 2 it is present for systems of finite size when 
the system size exceeds the mean free path but not the 
localization length. It is suggestive to conclude that dif- 
fusive motion is associated with a large spread in phase 
space and thus with the minimum observed in the inverse 
participation ratio. 

To substantiate this idea, we compare our results to 
energy level statistics. In the diffusive regime the en- 
ergy spacing distribution p{s) is close to the Wigner- 
Dyson distribution pw{s) = ^tts exp(— 7rs^/4). In con- 
trast, in the localized regime, p{s) approaches the Pois- 
sonian statistics pp{s) — exp(— s). To quantify the form 
of the distribution we evaluate 

lo Mpis) -Pw{s)) 

V = (6) 

Jo ds{pp{s)-pw{s)) 

which is particularly sensitive to level repulsion. Here, 
b = 0.4729 . . . refers to the first crossing point of the dis- 
tributions pwis) and pp{s). According to its definition, 
T] = 1 for a Poissonian spacing distribution and 77 = for 
a Wigner-Dyson spacing distribution. 

In the lower part of Fig. ^, ry is shown as a function 
of the disorder strength. For weak disorder, rj exceeds 
1 because of non-universal level statistics appearing in 
regular geometries in the ballistic regime. For strong 
disorder one finds Poissonian statistics as expected for 
the localized regime. At intermediate disorder strengths 
in d = 1, ry exhibits a minimum due to level repulsion, but 
no proper diffusive regime exists. On the other hand, for 
finite size systems in d = 2, an extended region is present, 
where the level statistics is close to that predicted by 
random matrix theory. This property is commonly used 
to identify a region of diffusive dynamics. For d = 3 the 
diffusive region survives even in the thermodynamic limit 
L — > 00. 

Fig. 1^ clearly shows that the decrease of the inverse 
participation ratio P is related to a plateau of rj at values 



close to zero. Using the spreading of the states in phase 
space as an indicator for diffusive behavior is therefore 
consistent with the results from level statistics. For d — 
3, a comparison between rj and the distribution of P even 
allows to identify the ballistic regime for weak disorder, 
where the inverse participation ratio essentially remains 
constant. 

In conclusion, we have shown that phase space proper- 
ties represent a new tool to describe disordered systems 
for arbitrary disorder strength. The inverse participa- 
tion ratio P in phase space provides information, which 
could not be obtained from the corresponding quantities 
Px in real space and Pfc in momentum space. It is only 
the first quantity which captures the appearance of a dif- 
fusive regime in two and higher dimensions. Given our 
demonstration that the physics in phase space is able to 
provide new insight into the dynamics of disordered sys- 
tems, these concepts likely will also make a valuable con- 
tribution towards understanding the challenging problem 
of the combined effects of interaction and disorder in few- 
body systems. 
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